WISC-TH-98-348 



THE SPECIFIC HEAT OF A FERMI GAS AT LOW 
TEMPERATURE: 
A closer look at the T InT behavior 



A. REBEIQ 

Department of Physics 
University of Wisconsin. Madison, WI 53706 

W.N.G. HITCHONQ 

Department of Electrical and Computer Engineering 
University of Wisconsin. Madison, WI 53706 



Abstract 

At finite temperature, a Fermi gas can have states that hold simul- 
taneously a particle and a hole with a finite probability. This gives 
rise to a new set of diagrams that are absent at zero temperature. 
The so called " anomalous" diagram is just one of the new diagrams. 
We have already studied the contribution of these new diagrams to 
the thermodynamic potential (Phys. Lett. A224,127 (1996)). Here 
we continue that work and calculate their effect on the specific heat. 
We will also calculate the finite temperature contribution of the ring 
diagrams. We conclude that the TlnT behavior due to exchange gets 
canceled by the new contribution of the new diagrams and screening 
is not essential to solve this anomaly. 
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In this letter we calculate the specific heat of an electron gas at low 
temperature and high density. This calculation is not just an improvement 
over the Gell-Mann calculation ||T|. It is fundamentally different. The Seitz 
formula used by Gell-Mann has been criticized in the literature and the rel- 
evance of screening has been questioned ||, |§. We would like to stress that 
we are not arguing that the divergence of the inverse of the density of states 
at the Fermi level is not due to the long range nature of the Coulomb field. 
Here instead we are addressing the original problem of the TlnT behavior 
of the specific heat that was pointed out by Bardeen || (see also ||) if ex- 
change is included. In fact our calculation shows that there is no need to 
invoke the screening of the Coulomb potential to find a linear function for 
the specific heat at low temperature. We stress that the solution to this 
problem has apparently not been clarified yet! Clearly finite temperature 
calculations show that the problem is related to the Fermi distribution. In 
fact all of these calculations || give a finite answer for the specific heat with 
exchange included. Unfortunately, Gell-Mann in his paper did not settle this 
point, even though he cited Bardeen's paper. This caused a confusion be- 
tween two different types of logarithmic behaviors. One is due to the long 
range behavior of the Coulomb field while the other is due to the statistic of 
the Fermi particles. In this letter we show that the specific heat maintains 
its observed linear behavior at low temperature without one having to worry 
about screening. However as we showed in a previous letter we have to 
include a new infinite set of diagrams. This new set does not contribute 
anything at zero temperature. The diagrams are merely due to the fact that 
at nonzero temperature there is a finite probability of finding a state that is 
occupied by a particle and a hole. 

We start as usual from the Hamiltonian of the system, 



Here, we use units such that h = m = e = 1. ty(x) is a two-component 
electron field, a and (3 are spin indices, i.e., a = 1 for spin up and a = — 1 
for spin down. Summation is implicit for repeated indices. The system is 
constrained by the condition 
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where N is the electron number operator which is constant. Because of this 
constraint, we prefer to work with the Hamiltonian 



where fi is a Lagrangian multiplier. This calculation we present is mainly an 
improvement on calculations of the exchange energy and ring diagrams made 
by Isihara and coworkers [fjj and others. The appearance of temperature de- 
pendent logarithmic terms in the internal energy found through various cal- 
culations are unacceptable. These terms will reappear again in the specific 
heat and will spoil the observed linear behaviour at low temperature. In the 
past these difficulties were not given close attention. They were either hidden 
in an effective mass or the expansion was not completely correct as we argue 
below. The history of this problem dates back to the '30s when it was first 
pointed out by Bardeen [[| . A preliminary treatment of this work appeared 
in J7|. It was argued that a correct approximation to the internal energy of 
a many-particle system up to exchange, i.e., a Hartree-Fock approximation, 
will automatically include summing an infinite set of diagrams (figure |T]). 
The first diagram in this set is known as the "anomalous" diagram ||. 
Therefore any approximations must be made at the level of two-particle 
Green's functions and not the single-particle wave functions as was done in 
most past treatments. This is exactly the reason for introducing two-point 
sources instead of the one-point source used in the usual effective potential 
method [|, @. We have seen that this enabled us to write an expression for 
the energy in terms of the true two-particle Green's functions of the system. 
Hence for the two-body Coulomb potential a Hartree-Fock approximation 
is equivalent to writing the two-particle Green's function in terms of the 
one-particle Green's function as follows: 



This approximation is equivalent, as it is well known, to the usual Hartree- 
Fock approximation of the energy to order e 2 at zero temperature. At low 
temperature there is another contribution to order e 2 that comes from sum- 
ming an infinite set of diagrams that were never treated before. The ring 
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G(xi,x 2 ; x 3 , x 4 ) = G(xi, x 3 )G(x 2 , x 4 ) — G(xi, X4)G(x 2 , x 3 ) 
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diagrams treated in [|I] are of order e 4 and therefore the new set can not be 
ignored in any finite temperature treatment. The effective action method 
can be easily extended to the non-zero temperature case. In the presence of 
external sources, the grand partition function Z of the system is given in this 
case by the following expression: 

Z = exp(-/30) = Tr exp (-/3 (H - jxN) + j(p + -<pB<p + ^ f Q^ ) . (5) 

The trace is over space indices, internal indices and "imaginary time" . The 
imaginary time interval is [0,/3\. The thermodynamic potential Q should be 
considered as a functional of the sources j(x, r), B(x, r; x', r') and Q(x, r; x r , r') 
The other parameters in the system are the variables p, T and V. Now we 
define new variables by the following functional relations: 



JJ = Vc, (6) 
— = A + (7) 

and 

The above equations are evaluated at j — B — Q — 0. To get the free energy 
F[p, <p c , A] we make the following triple Legendre transformation 

F[p, <p B ,A]-nN = B, Q] - jip c - hp c B<p c -pQ-AB (9) 

To calculate F perturbatively, we expand the exponent around the classical 
Hartree potential in a neutral background. The coefficients of this expansion 
are expressible in terms of Feynman diagrams. We keep only the first two 
terms of the expansion. The equations of motion of the variables p, ip c and A 
are obtained by the standard way. The variables ip c and A can be expressed as 
a series expansion in e 2 with coefficients dependent only on p. Hence within 
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the above approximation, the expression for the thermodynamic potential 
has the following form: 



1. 
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A = Vl + Vl ex + -Tr {In 5(x - y) + e 2 ^(x - y) - e 2 $l(x - y)} 



— Tr{ln 5(x — y) — e ^s(x — y) + e ^s(x — y)} (10) 
where the functions $l(x — y) and Q(x — y) are given by 

$l(x -y) = [ dr I d 3 zA(y - z)p(x - y)p(z - x) (11) 

and 
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S(x — y) = I dr I d 3 zA(x — z)p{z — x)p{y — z). (12) 
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The term Q is the free contribution, Q ex is the usual first order exchange 
term, the third term represents the usual ring diagrams and the last term 
represents a new series of diagrams shown in figure 0. The function A(x — y) 
is the bare Coulomb potential and p(x — y) is the free one-particle Green's 
function at finite temperature. The signs in front of the two Tr ln-terms 
are determined by integrating over Bose fields and Fermi fields, respectively. 
We have shown that the new diagrams provide a T In T term with the right 
weight to cancel the TlnT term that appears in the exchange term. This term 
is beyond the reach of all previous treatments including those of Gell-Mann 
and Isihara. Hence the unphysical logarithmic behavior of the exchange 
energy gets cancelled by a consistent perturbative treatment of the many- 
body problem. The other important point we learned from this calculation 
is that this logarithmic term is not an immediate byproduct of the singular 
behavior of the Coulomb potential since these new diagrams are not singular 
like the ring diagrams. The new contribution at high density is given by 

' ( ' V ;kUC + hn V ). (13) 
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The ring diagrams contribute a temperature dependent term of order e 4 
as we will show next. The ring diagrams are a manifestation of the long 
range behavior of the Coulomb field. In fact, in our calculation, the ring 
diagrams are obtained by integrating over the Hartree field while the new 
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terms are obtained by integrating over the Fermi fields. We do not think 
that this point has been recognized in the literature. Obviously, these new 
diagrams have a nonzero contribution to the specific heat and must not be 
left out in any accurate calculation. Their contribution appears when we try 
to calculate fi in terms of the Fermi energy of the free electron gas. The fact 
that these new terms have an opposite contribution to the exchange energy 
can be easily deduced if we refer back to [0. In fact, our method shows that 
for a many-body problem with a two-body potential, isolating the exchange 
energy is really an artificial effect. Reference [|K| shows that the exchange 
energy term really comes from the ring diagrams by integrating over the 
Hartree potential. However, the new terms came by integrating the Fermi 
field degrees of freedom which have a statistical origin. 

Now we continue evaluating the contribution of the ring diagrams at finite 
temperature. This contribution is given by 

R = ^Tr{ln[5{x - y) + e 2 K(x - y)\ - e 2 K(x - y)}. (14) 

In momentum representation, this is given by 

V n=co 4-tt r°° 
R =^Htt^ p 2 ^[ln{l + e 2 ^,c„,)}-e 2 3?(p,u; n )] (15) 

where p represents momentum and u n = ( 2n + 1 ) 7r ; n — 0, ±1,±2, .... The 
function 3ft(p, u n ) is given by: 

*<» ^ - ^J ^wfv%* m ~ np+ ql) - (16) 

The function f(q) is the Fermi-Dirac function, 

1 + exp[^(g 2 /2 - n)\ 
To evaluate R we break up 3?(p, cu n ) in two pieces. We write 

n(jp,u n ) = F l (p,u n ) + F 2 (p,u n ) (18) 
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such that 

= ^ J d \ E *!l%rU m (19) 

where Eg- = y. The other function F 2 (p,uJn) is similarly defined but with 
f(q) replaced by f{p + q). 

We start by evaluating the function Fi(p, uj n ). Taking p along the z-axis and 
integrating the angular coordinates, we get 

F( 1 [°° , 1 ^l+(pg + p 2 /2) 2 1 

r l{ p,u n} (2yr)V y Q qaq m^ + {pg _ p2/2)2 j 1 + exp ^ q 2 /2 _ ^ " 

(20) 

Now we rewrite this in the following form: 

Fi(p,w„) = / grig In 2 _ 7o\2 

{2n) 2 p- i Jo u)^ + (pg - p z /2y 

1_ ^ + (pg + p 2 /2) 2 



+ 



(27r)Vy * * L w 2 + (pg-p 2 /2) 2J 1 + exp{-/3(g 2 /2-/i)] 
(2n)2p3Jk F q q [ cu 2 + ( W -p 2 /2)2 J l+exp{/3(g 2 /2- /U )} • 1 J 



In the limit when /3 — > oo, the first integral that we isolated becomes simply 
part of the zero temperature contribution to R. We denote this part by 
F lfi (p,u n ). Next we set z = f3(q 2 /2 — /i) and integrate over z instead of q. 
Since we are interested in the low temperature limit, we get to leading order 
in 1/(3 that 



7f 1 roo 

Fi{p,u n ) = F lfi (p,u n ) + . , - / — {In 

(2nYv 6 3 Jo 1 + exDz 



7T 1 too dz Ujl + {p 2 /2+pk F {l + ^)f 

(2tt) V /3 7o 1 + exp ^ ! n[ cu 2 + (p 2 /2 - pfc F (l + -M) sJ 



cu 2 + (p 2 /2+pMl-p 2 
L 2 + (P 2 /2-pA; F (l-f)) i 



- M ,o I I 2 / o 1 , ? J ) • ( 22 ) 



Again this last integral over z can be separated into four integrals that can 
all be evaluated in the same way by elementary methods. To leading order 
in (3, we find 



1[P ' Un) " ^ + (2vr)2pA; F 12/3 2 ^ + (p 2 /2 + pM 2 + ^ + (p2/2 _ pkp y > 



The function F 2 (p, oj n ) is easily seen to be equal to Fi(p,u n ), hence to this 
approximation the function 9ft (p, u; n ) is equal to 

£(p, co n ) = 2F 1)0 (p, cu n ) + 2Fi(p, cu n ) . (24) 

Putting his back into i?, and expanding the logarithmic function in powers 
of i, we get 

V An n=OD r°° 
R = ^7^ E / P 2 ^ln{l + 2e 2 F 1 ,o(p,cu„)}-2e 2 F 1 ,ob,^) 

2/? (2tt) 3 n fr' 00 /5 2 Jo 1 + 2e 2 F lfi (p, w„) 

where we have defined cj„) to be 

x 1 r 9 + ^ p/2-k F , 

24p/c F + Qtr/2 + pk F y uj 2 + \p /2 - pkf) 

The first integral in i? will be denoted by i?o an d the second one by Rp. So 
far our approximation works for any density. Now we notice that F lj0 (p, u> n ) 
is proportional to the inverse density. Hence for high densities the second 
integral can be expanded in powers of inverse densities. To lowest order, the 
second part can be easily evaluated. Letting x = and y n = pf , we have 

e 4 V " r°°dx ( x/2 + 1 x/2-1 \ 

~ To7oZ\4u3~m ~n ^ L \ „,2 i „2C„/o i ^2 ^ 



12(2tt) 4 A;3/3 2 13 n ±?Jo x 2 \y 2 n + x 2 {x/2 + l) 2 y 2 n + x 2 {x/2 - I) 2 ) 

r x(x/2-l). r x(x/2-l), 
{x — y n arctanj } — y n arctanj } 

+ 1/2 h /" + X2{X/2 + 1)2 1 + -^(4, 2 - x*) \n\ vl + X2{X/2 + 1)2 ] j (27) 
+ i/Z[n[ y 2 +x 2 ( x /2-l) 2l + 8x 2 ^ Vn X )ln[ y 2 +x 2 (x/2-l) 2n - W 



In the limit that (3 — > oo, the summation over n can be changed to an 
integral over y with measure dy = -J^-. This integral has been evaluated 

numerically. It turns out that integrating over the variable - instead of y 
makes the integration more suitable from the numerical point of view. The 
numerical value of the integral 



I = — dzK(x,z)P(x,z) (28) 
Jo x z Jo 

is approximately equal to —22.94. In the above the functions K(x, z) and 
P(x, z) are given by 



K( \- 1+x / 2 l-x/2 
{X ' Z) ~ z 2 + (1 + x/2) 2 ~ z 2 + (1 - x/2) 2 { } 



and 



„, , , fl + x/2\ (l-x/2" 
P(x, z) — 1 — ^arctan — ^arctan 



z 

+ J_ (i _ x 2 / 4 + Z A ln[ f + / (1 + /2) J . (30) 
2x^ 1 ' l z 2 + (l-x/2) 21 K 1 

Therefore the finite temperature contribution of the ring diagrams to leading 
order in 4j- and for high density is given by 

R = iV(0.0622 lnr s - 0.142) + 0.0181 N(T)fi( aTs ^ ) 2 (31) 

Here N is the total number of particles, which is a constant. All the other 
variables are defined as follows: 



/ 4 me 2 h 2 

« = (tH 173 ,^) = — 2— , a 
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Note that here we did not resort to any artificial cutoff of the momentum 
in our calculation as was done by Gell-Mann |1J or Isihara 0. The first 
author assumed T = 0, while the second author and coworkers carried out 
the integral up to p < 2kp- In fact momenta larger than 2kp contribute 
about 5 percent to the answer. 

Our calculation of the specific heat starts by finding an expression for \i in 
terms of the Fermi energy at zero temperature. In other words, we expand \i 
in powers of Jy. We need to do this before any differentiation with respect to 
temperature, to get the free energy function which is by definition a function 
of T, V, and N. Here we will need to carry our expansion only to first order 
in -k. This will at least give us the first nontrivial contribution beyond the 
ideal one. Hence in this case it will be enough to use for /i the equivalent 
free electron expression, that is we take, 



71 



(33) 



This implies that 



N(T) = N(l - -5) . (34) 

However for r s (T), it is enough to take it equal to the corresponding zero 
temperature value r s . The calculations that lead to the determination of the 
temperature dependent terms in the free energy are straightforward. The 
ring diagrams will contribute a term Rt, which in Rydbergs is given by 



N r , vr 2 w , 22.90 . . , 

R T = — (1 0.06221nr s - 0.142 + — — - . 35 

^ LV 8 A 1 64(2tt) 2J v ; 

Similarly, the "anomalous" diagrams contribute a term of the form 



3vriV _ 1 
2r\l(xr s 



A t = -tt^— (C + T^lnr/o)- (36) 
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The constant ( is approximately equal to 8.0. Besides these, we have the 
usual contribution from a kinetic part and an exchange part. They are, 
respectively, given by 



Kr = 2 -^K (37) 



and 



The constant a has different values in the literature. For example Isihara 
and coworkers give it a value of approximately — ~. Our calculation gives it 
a value of almost +2.0. Both values give similar answers to the specific heat. 
It will be seen below that the contribution from the new diagrams is much 
more important and hence any discrepancies in the values of a will hardly 
matter. 

Remembering that F(T, V, N) = Q + fiN, we can easily find the contribution 
of all of the above terms to the temperature dependent term, Ft, of the free 
energy F using standard techniques ||11|| . Finally, we can now get an estimate 
for the specific heat , 

d 2 F 

Cv = -T-^ (39) 



such that 

F(T, V, N) = F (T,V,N) + : 0.077 In r s + 0.184 (40) 

fx \ r s i 



and Fq is the free energy of a free electron gas. Therefore, if Cy is the specific 
heat for the ideal Fermi gas, we get for the specific heat at low temperatures 
and high densities the expression, 

C 

= 1.00 + 0.185r s -r s 2 (0.0017- 0.00071nr s ) . (41) 

Cy 
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Some comments on this expression are called for. First, without the new 
effect of the anomalous diagrams, the second term would have been negative 
besides the InT dependency. Second, for sodium, which has r s = 3.93, the 
ratio Cy/Cy = 1.71 compared to the experimental value of 1.26. This is not 
a bad outcome considering that our calculation is strictly valid for only small 
r s . Gell-Mann's computation gives 0.48 as an answer in this case. Going 
back to Eq.(22), it is easily seen that had we kept the next term in our 
approximation, we could easily improve our answer. Besides the numerical 
answer, it is important to realize that screening has little effect on the specific 
heat and need not be included to get a finite and reasonable answer for the 
specific heat. Exchange energy is essentially of statistical origin. At nonzero 
temperature a point in phase space can simultaneously hold a particle and 
a hole and hence the statistics of holes becomes important as well and it 
is manifested in the appearance of this new type of diagrams. The name 
"anomalous" is therefore misleading. In fact, they should be called particle- 
hole exchange diagrams. This name shows that this kind of exchange is not 
allowed at zero temperature. 

We end this letter by stressing the fact this expansion is not an expansion 
in e 2 , and that what this method considers exchange is not the usual meaning 
of the word [TIJ. The usual first order exchange, the ring diagrams and 
the new particle-hole exchange diagrams are all part of the first correction 
term to the energy of an interacting Fermi gas. The second order exchange 
diagram and the ring diagrams with exchange appear in the next term of the 
expansion. The third term in this expansion contains among other things 
the ladder diagrams. These diagrams are important for metal densities. We 
hope to include all these effects in another communication. 

Finally, we would like to mention that this expansion with two-point 
functions was also recently used in the study of the anharmonic oscillator at 
low temperature ||12|| . Here it was shown that effective actions with two-point 
functions give better results than effective actions with one-point functions. 
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Figure 1: A new infinite set of "anomalous" diagrams. In this approxima- 
tion, the solid line corresponds to the Fermi propagator and the dashed line 
corresponds to the Coulomb field. 
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